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Abstract: We generalize (linearized) relativistic hydrodynamics by including all order gradient 
expansion of the energy momentum tensor, parametrized by four momenta-dependend transport 
coefficients, one of which is the usual shear viscosity. We then apply the AdS/CFT duality for 
J\f = 4 SUSY in order to compute the retarded correlators of the energy- momentum tensor. 
From these correlators we determine a large set of transport coefficients of third- and fourth-order 
hydrodynamics. We find that higher order terms have a tendency to reduce the effect of viscosity. 
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1. Introduction 

Relativistic Heavy Ion Collider (RHIC) produces hadronic matter with temperatures ranging between 
the initial Tj ~ 2Tc and the final (or freezeout) value Tf ~ Tc/2 ^, where Tc ~ 170MeV is the 
QCD critical temperature. It has been shown fs], ^ that the flows (radial, elliptic) associated with 
the plasma expansion are well and consistently described by near-ideal relativistic hydrodynamics, 
with freezeout implemented via hadronic cascades. Since elliptic flow is dominated by the early times 
of the fireball expansion, when T > Tc and matter is in the so called Quark-Gluon Plasma (QGP) 
phase Q, this lead to the conclusion ^ ^] that QGP is a "perfect liquid", presumably because it 
is actually in a new - strongly coupled - regime of QCD. Those views were discussed in detail and 
eventually accepted in the 2004 "white papers" of all four experimental collaborations Q. The full 
understanding of QCD dynamics at strong coupling, even in the deconfined phase, remains a challenge 
and one usually appeals to either lattice simulations or phenomenological models. While the lattice 
is considered a reliable source for QCD thermodynamics, it usually fails to provide accurate data on 
transport coefficients. Thus, in order to understand transport properties of QCD, at the moment we 
have to appeal to various microscopic models (for recent reviews on the strongly coupled QGP see 
e.g.|l,|ll). 
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Starting from 2004 RHIC experiments have discovered and studied phenomena known as the 
"cone" and the "ridge" , associated with the propagation of the energy deposited by quenched hard 



jets (for the description of the phenomena and recent data see Refs. [12, ^ and references 
therein). The former was associated with conical hydrodynamic flows induced by fast particles 
propagating through the medium ||l^. This has in turn initiated studies, within AdS/CFT, of such 



processes induced by a steadily moving heavy quark, see Although these studies 

did not appeal to any hydro description, their results were found to be in very good agreement with 
(even unimproved) hydrodynamics. 

Another structure, known as "soft ridge" has been observed in two-particle correlations: its 
origin is attributed to initial state fluctuations in the colliding nuclei. For experimental data and 
phenomenological discussion see the talks at the BNL dedicated workshop |2^ . Although in this paper 
we will not discuss the phenomenology of those objects, we nevertheless stress that they provide 
the strongest motivation for a detailed study of small perturbations on top of (hydrodynamically 
expanding) matter. 

One obvious step in understanding these perturbations is to study them using linearized hy- 
drodynamics, assuming their amplitude to be small. On the other hand, these objects start their 
evolution at much smaller scales compared to nuclear (or fireball) radii. For example initial state 
fluctuations are believed to be given by the "saturation scale" 1/Qs, which is only 0.2 fm, about 30 
times smaller than the fireball as a whole. Therefore, the evolution of small perturbations includes 
much larger spatial gradients, and in order to treat them better one would naturally try to improve 
the accuracy of hydrodynamics, including higher order derivative terms. Since the latter appear with 
many new transport coefficients, the usual phenomenological approach which derives viscosity from 
the data would hardly be possible. Instead, some self-consistent approach is needed to calculate as 
many of them as necessary. 

Such model of choice for the present study is = 4 SUSY at large Nc- Via the celebrated 



AdS/CFT correspondence |21| this gauge theory at strong coupling admits a dual description in 
terms of weakly coupled gravity in AdSs S5 space. The finite temperature version of this field 
theory is dual to the AdS-Schwarzschild black hole (brane). The laws for Schwarzschild black hole 



thermodynamics imply that the entropy density is proportional to the area of the horizon |23] ^. The 



equilibrium pressure is P = vr^ A^^T^/S, while the energy density is e = 3P due to the conformal 
symmetry of the microscopic theory. 



Refs. 1 24 1 pioneered the study of transport coefficients via dual description. For a static plasma 
and in the limit of large 't Hooft coupling A 3> 1, the ratio of shear viscosity to entropy is independent 
of the coupling and is in fact remarkably small 

s Air 

Furthermore, Refs. [24| conjectured that this value for the ratio is a universal lower bound, valid for 



all physical systems in nature. While AdS/CFT leads Eq. ( |l.lD , it does not provide any explanation 



^For theories involving higher order curvature corrections to the Einsteinian gravity, the relation between area of 
the horizon and entropy is invahd p3]. 
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of this result from the gauge theory side. So far, no microscopic mechanism for low viscosity has 
been established for QCD, though a promising proposal can be found in [p5| . 

It is rather difficult to extract viscosity from experimental data precisely, because it is small 
and its effects are 0(10%) or so, comparable to other uncertainties. The phenomenological studies 
of RHIC data (such as in Ref |2^) typically focus on the elliptic flow dependence on centrality or 
transverse momentum, V2{b,pt)- The optimal ijq/s for these fits occurs at a value of the order of the 
suggested minimum, although deviations from it by a factor two or so are still possible. Another 
argument to support very low viscosity comes from discussions of the overall entropy production, such 
as in Ref. [^. Those works suggested that there is a tension between the total entropy (measured 
by the observed multiplicity of produced hadrons) and the very short thermalization time (initial 
time for hydro evolution), unless the viscosity over entropy ratio is pushed down, maybe even below 
the bound. The third (more indirect) argument for low viscosity is the survival till freezeout of 
the "cones" and "ridges", suggesting smallness of dissipative effects. Therefore all these approaches 
indicate a very small viscosity value. 

The hydrodynamic representation of the energy-momentum tensor is 

{Tf"") = {e + P) + P g^"" + U^" (1.2) 

where the average is taken over the thermal bath. While at some microscopic scale I the system 
is assumed to be locally at thermal equilibrium, at some macroscopic scale L ^> I the local fluid 
velocity field u is a function of space-time coordinates. The "tensor of dissipations" 11'^'^, added to 
the ideal- fluid part, represents all the deviations from the equilibrium state induced by such a flow 
field. In the long wavelength limit L ^ I, these fluctuations can be expanded in terms of gradients 
of the velocity field, or in powers of l/L. The first order Navier-Stocks (NS) hydrodynamics retains 
only the first gradient ^ 

n^'^ ~ ryoV^u^ (1.3) 

In this work we will discuss higher order gradients, which will provide certain corrections to the first 
order viscosity term when gradients grow. These corrections are relevant for smaller size objects in 
the plasma or to earlier time of hydro evolution. 

The high order gradient expansion generically includes two types of terms: (i) non-linear terms in 
the velocity field (like (Vn)^) and (ii) linear terms with multiple gradient operators acting on a single 
velocity field (like VVu). These two types of terms are controlled by two different parameters. The 
non-linearities are important when the field amplitude is large. However, even for small amplitude 
waves, one can get large contributions from the linear terms when the momenta associated with the 
wave are large. 

Recently, second order hydrodynamics (next-to-NS) attracted significant attention [||, |2|. The 
main reason is that NS hydrodynamics is known to have causality problems. The accausal effects 
create numerical instabilities when solving hydrodynamic equations. The problem originates from the 
fact that NS equations imply instantaneous response to any perturbation introduced in the system. 

^Throughout this paper we will be considering conformal theory only, for which there is only shear viscosity since 
the bulk viscosity is zero. 
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In order to circumvent this problem, one may introduce a relaxation time. It explicitly appears as a 
new transport coefficient when the gradient expansion is extended to second order: 



~ r/o [1 - r V) ] V^n'^ ^ ?/o [1 + ir to] (/t^ u") . (1.4) 

In fact, in order to restore causality it is not sufficient to include second gradients only: all order 
gradients need to be resummed. A very popular resummation scheme is due to Israel and Stewart 
(IS) It essentially generalizes viscosity to an cj-dependent but /c-independent (complex) function 

""H^^. (1.5) 

Eq. (^]^) can be viewed as a Pade-like resummation of ( |1.4| ). The relaxation time r provides a scale 
for exponential relaxation. The position of the pole below the real axis and the "good" falling off 
asymptotic behavior make the model causal. In other words, as a function of complex ui, the viscosity 
is analytic in the upper half plane. In coordinate space this simple pole corresponds to a memory 
function with exponential falloff: 

n'^^(x,t) ~ ^ r dt'e-(*-*')/" V^n'^(a;,tO. (1-6) 
Jo 

In this paper we will be studying all order gradient expansion in the linear approximation. Instead of 
introducing new transport coefficients at each new order, we will be thinking of viscosity and other 
transport coefficients as frequency and momentum dependent functions. We will be working in the 
framework of AA = 4 SUSY. For the rest of this paper we set all dimensionfull units to be related with 
the temperature, 2itT = 1. Among our results, we will show that the IS resummation, although 
well-known and used, is still simplistic model for high order gradient terms, which is even qualitatively 
inconsistent with AdS/CFT results. Not only it misses important non-linear terms already at second 



order |2£, but (as we will show below) it is also incorrect in the linear approximation starting 



from the third order. 

More generally, we will find that higher order terms do have a tendency to cancel (or reduce) 



the effect of NS viscosity. In particular, in our earlier paper |31|, we argued that the extremely low 



viscosity suggested by Refs. ||2g, ^ may essentially be some "effective viscosity", which includes 



these high order gradient terms. The real systems probed in RHIC collisions have finite gradients and 



the inclusion of their effects may demand going beyond NS approximation. In |31] we attempted to 
extract a momentum-dependent viscosity from the imaginary part of the sound dispersion curve. Our 
main observation was that the effective viscosity as probed at finite momenta turns out to be smaller 



compared to the value at the origin. Motivated by [^], we discussed in |31] the implications of a 
momentum-dependent viscosity on the entropy production for Bjorken expansion |^^. We discovered 
that the inclusion of momentum-dependence made it possible to push the hydrodynamic description 
a bit further into earlier times of the collisions, with the entropy production due to viscous hydro 
stabilized at around 20% of the total entropy produced in the collision. The conclusion is that the 
account for a momentum-dependent viscosity reduces the sensitivity to thermalization time. Now, 
with the result reported below, our previous approach [^] based on the sound dispersion curve looks 
rather naive (for a much more elaborated study of hydrodynamic theory as an effective theory for the 
lowest modes see Ref. p3|). In general, we will see that the sound dispersion curve does not contain 
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enough information to define the "generahzed viscosity" function. Nevertheless, we qualitatively 
captured at least the right trend: full second order hydro with all non-linear terms included has the 
same trend towards reducing the entropy production 34 ] . 

Our goal in this paper is to put the idea of a momentum-dependent viscosity on a more solid 
ground compared to our naive treatment in |3l|. In the present analysis we will be focusing on the 
retarded correlators of the stress tensor. The correlators contain information not only about the 
positions of the poles but also about their residues. The complete information on the correlators is 
equivalent to the knowledge of the energy momentum tensor in the linearized approximation. 

In a conformal theory in four dimensions, there are only three independent correlators of the 
energy-momentum tensor. These are correlators in the sound (G^), shear (G^), and scalar {G^) 
channels. AdS/CFT correspondence provides a tool to compute these correlators by solving certain 



linearized gravity equations in the background of the AdS-Schwarzschild black hole p5| , pq , 37|. These 
equations essentially describe graviton's propagation from the AdS boundary, where the field theory 
is defined, to the horizon of the black hole. Absorptive boundary conditions are imposed there. 
Dissipation takes place at the horizon while there is no dissipation in the bulk of AdS. However, 
the bulk curvature acts as a non-linear medium, which provides a source for complicated dispersion. 
It is this dispersion which, by means of the duality, is mapped into momenta-dependent transport 
coefficients. 

Our strategy is to first write a most generic hydro-like representation of the energy momentum 
tensor T^*^, in terms of the fluid velocity fleld u. We find that, generically, there are four structures 
(or operators involving derivatives of u or the metric g) which can occur in T'^'^ and are consistent 
with all symmetries. Each structure enters with a coefficient which is momentum-dependent. These 
are the generalized transport coefficient we are looking for. One of them is associated with the 
shear viscosity, while the remaining three encode responses of the system to external (4d) gravity 
perturbations. We call them gravitational susceptibilities of the fiuid (GSF). The operators which 
are multiplied by the GSFs involve the Weyl tensor of the metric and vanish in the fiat Minkowski 
space. 

We the proceed by using this hydro-like representation of T^'^ in order to compute its correlators 
in the three channels introduced above. We then attempt to determine the momentum-dependent 
transport coefficients from the matching to the functions G^ , G^ , and G^ computed directly from 
the bulk gravity side. 

Our program runs into a problem, which we were not able to resolve completely: there are in 
fact four independent transport functions to be extracted from three equations. Despite the fact that 
we could not determine the entire functions, we were able to get them to quite high order in the 
perturbative expansion at small momenta. In particular, we found the shear viscosity function to 
fifth order in the gradient expansion. This involves several new transport coefficients, most of which 
are obtained numerically. 

The conceptual problem mentioned above, prevented us from computing shear viscosity in the 
whole kinematic region of arbitrary frequency and momentum. Instead, we build a model similar 
to IS which utilizes the information about the new transport coefficients and preserves the causality 
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condition. We propose this model for phenomenological studies of hydrodynamics at RHIC, but any 
application of this model is left beyond the scope of this paper. 

The paper is organized in the following way. In Section 2 we present the general setup for 
computing the retarded correlators from the bulk gravity and from the generalized hydro on the 
boundary. Section 3 presents some results. A phenomenological model for generalized viscosity is 
proposed in Section 4. Our conclusions are summarized in Section 5. Two Appendices supplementing 
Section 2 provide details of some analytical computations. 

2. Generalities 

The retarded correlators of two energy-momentum tensors are defined as follows 



Here the average is over the equilibrated thermal bath. For conformally invariant plasma with 
traceless T^^ , there are only three independent correlators = G^^^^ (tensor), G^ = G*^'*^' (shear), 
and G^ = G*^*^ (sound) with the vector k pointing in the z-direction. All other correlators are related 
to these three either by rotational symmetry or by the equations of motion. 

2.1 Life in the bulk: Retarded correlators from gravity 

In this subsection we closely follow the setup and results of Ref. |^]. From the bulk gravity 
side, in order to compute the retarded correlators at non-zero temperature one has to solve certain 
wave equations (one for each symmetry channel). These equations describe propagation of the 
corresponding metric perturbations (gravitons) in the AdS-Schwartschild BH background of the dual 
description. The differential equations are of the form 



where the coefficients Pair), Qai''') depend on the frequency a; and momentum k, and a = T, D, S 
labels the three symmetry channels. The coefficient functions are given by the following expressions. 

• The scalar channel 




(2.1) 



— ^Za(r) +Pa{r) — Za{r) +qa{r)Za{r) = 



(2.2) 





(2.3) 



where / = 1 — r 



The function / is inherited from the AdS-BH metric. 



• The shear channel 



{io^ - k^f)f + 2r^io^ 
r/(A;2/-u;2) 




f 



(2.4) 



• The sound channel 



ps{r) 



3a;2(l + r^) + k^{2p - 3r^ - 3) 
r/(3w2 + ^2(^2 _ 3)) 



qs{r) 



3a;4 + A;4(3-4r2+r4) + k'^{APuj^-6uJ^-Ar^ f) 
r/2(3cj2 + A;2(r2-3)) 



(2.5) 
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The fifth dimension coordinate r ranges from to 1, where r = corresponds to the boundary of the 
asymptotically AdS space, and r = 1 corresponds to the event horizon of the background metric. 

The information about the retarded correlation functions is encoded in the solutions to Eq. ( |2.2| ), 
which satisfy the incoming wave condition at the horizon Za{r — > 1) ~ exp[—iuj/2]. At r = the 
solution can be written as a linear combination of two independent local solutions, 

Za{r) = AaZiir) + BaZi\r), (2.6) 

Here Z^ is irregular in the origin while Z^^ is a regular solution. 

The prescription to compute the correlators G follows from the Minkowski formulation of the 
AdS/CFT correspondence and amounts to computing the ratio between the two coefficients in the 
expansion (l2.6|) 

G%u:,k) = -SP^f^ (2.7) 

For the three symmetry channels the correlators G are related to G, 

1 1 A-2 

Qxyxy = _qT. Qtxtx ^ _ qD ^3.8) 

2 2 u)'^ — k^ 

For the sound channel the relation is a bit more involved and includes a contribution from contact 
terms |p9| 

k^ ~a 1 29k^ - SOk^uj^ + 9uj^] 
■G^ + — 



r-itttt _ 1 
^ - 2 



3 (w2 _ ^2)2 - 12 (A;2 



(2.9) 



Eq. (^]^) has real coefficients which are even functions of frequency. In other words this equation 
propagates waves without any dissipation. The dissipation (time irreversal) effects are introduced by 
the boundary conditions at the horizon. However, the AdS-BH metric acts as a non-linear medium 
for the propagating graviton. The non-linear dependences on frequency and momenta which appear 
in ( p.2| ) are to be mapped onto highly non-trivial momenta dependence of the transport coefficient 
functions. 

Let us make a technical remark on numerical solution. The equations for the shear and sound 
channels have singular points inside the bulk r = [0, 1]. For the shear channel it appears for cu < k 
at ro = -s/l — io'^/k'^ and for the sound it is at tq = y^S (1 — uj'^/k'^) (condition that ro is inside the 
bulk). It would be interesting to understand if these points have any special physical role. To ensure 
that there is no instability caused by these singularities, we split our numerical solution into two 
intervals [0,ro] and [ro, 1] and matched the solutions at the singular points using analytic solutions 
in the vicinity of tq. 

The correlators computed from the gravity side agree with the field theory correlators up to a 
constant [35|. In particular, in the sound channel the relation between the correlators G**** and G*^*^ 



IS 

,2 fr^tttt , \ _ T 2 tr^tztz 



(C"" + e) = k' (G*^*^ + P) (2.10) 

The analytical expansion for the correlators at small momenta can be found in Appendix A. For 
the shear and sound channels we show some numerical results alongside the corresponding curves for 
the NS and IS hydrodynamics on Figs. || and ^. 
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Figure 1: Shear channel: top k = 0.4; bottom cu = 0.4. Solid line corresponds to the AdS/CFT correlator. 
Short dashes display the NS hydrodynamics while long dashes show the IS hydrodynamics. 

2.2 Life on the boundary: All order hydrodynamics 

The thermal field theory on the 4d boundary is defined by means of the generating functional 

Z[h] = I L»</.e^«[^l + Sd'-h.^T^" (2.11) 



where (j) collectively denotes all fields of the theory, Sq is the fiat metric action and h^u is an external 
perturbation of the Minkowski space. 

The expectation value of the energy-momentum tensor at non-vanishing external field h^y is 

(51nZ 



6h 



(2.12) 



Within the linear response theory we keep terms linear in h only. The correlators G'^^^'^ differ from 
the retarded correlators G"^^'^ by constant contact terms p5[ . 

We use eq. ( 2.12| ) to define hydrodynamic variables. Here {T^^'^)^^^ corresponds to the thermal 
equilibrium. The equilibrium energy density is 

60 = {T'Xr' (2.13) 

The external perturbation h^y shifts the theory from its thermal equilibrium. The out-of- equilibrium 
energy density is 

6 ^ (tOO)^^^ = eo + /i„^G"/3oo (2.14) 
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Figure 2: Sound channel: top k = 0.4; bottom u = 0.4. Solid line corresponds to the AdS/CFT correlator. 
Short dashes display the NS hydrodynamics while long dashes show the IS hydrodynamics. 



We can also define fluid's three-velocity 

(60 + Po)v' ^ {T'Yci = VG-'30* (2.15) 



and fluid's 4-velocity n'^ = (Vl + v'^, v) satisfying = —1. 



The action (2.11) is required to be invariant under the local Weyl transformation (see extensive 
discussion in Ref. ||29|| ) 

g,. - e-2^(-'*)5/... (2.16) 
The invariance of the action implies that T^^'^ and the velocity field u transform homogeneously 

T^^ g6n(x,t) ji/xi/ . _^ ^ (2.17) 



Weyl tensor 



In our construction below we will be imposing the Weyl invariance. To this goal we will employ the 

AA^A A A A-'-A A 

^liua ~ ^ ij.ua ~ 2 ^^i^ -^^ict ~ da ^liu ~ 9ijlp ^ a ~^ 9fia R y) ~\~ '^R {Qu 9fJ.a ~ da 9fju), 

which is constructed to be invariant under this transformation. Here R'^^^a^ R^ R stand for the 
Riemann, Ricci tensors and the scalar curvature. 
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The hydro representation of the energy-momentum tensor 

{Tf"')^i = (e + P)u^u^ + Pg^"" + H^f""^ . (2.18) 

Here for any tensor li^^ we define its traceless and symmetric component (following the notations of 
Ref |2|) 

^la,A 1 . . U.B ^ ^^^^^ _ _ ^^^^ 



n(M^> = i A'^^ (n«^ + n^„) - \ a^^ a"^ (2.19) 

with the projector 



^i^u = gi^^ + ui'u'' (2.20) 
which is commonly introduced to ensure transversity of the tensor of dissipations 11: 



u^W = 0. (2.21) 



This transversity is equivalent to the condition of no dissipation in the fluid's rest frame. The 
tracelessness of T^'^ implies that n^'^'^^ is also traceless ^ . The metric g is the full metric, namely the 
Minkowski metric perturbed by h. 

The energy-momentum conservation leads to equations of motion for the fluid: 

V^iTf"") = 0. (2.22) 

Here stands for covariant derivative with respect to the metric g. 

The tensor H^'^ is considered to have all order gradient expansion. Within the linearized ap- 
proximation discussed above, and constrained by the Lorentz and Weyl symmetries, there are four 
independent structures (operators) one can write down ^ 

= -2t]V''u'' + 2KUaUf3 C^''''^ + p {Ua V /3 + U/3 V„) C"^"^ + ^ V„ V/3 C^'"''^ . (2.23) 

By representing H^'^ in the form ( |2.23| ) we essentially postulate a constitutive relation between (T*-') 
and v^. This structure implies that the fluid can be perturbed either by inducing some velocity 
perturbation or by shaking the metric. These perturbations are not fully independent and can be 
related by the equations of motion: the gravity perturbations create perturbations of velocity (see 
Appendix B). 

Each of the four transport coefficient functions rj, k, p and ^ are considered to be functions of 
the Lorentz scalar operators and (u V) 

r? = r?[V2, {uV)] ■ K = k[V^, (uV)] ; p = p[V^ (nV)] ; ^ = C[V^ (uV)] ; (2.24) 

In momentum space representation (adequate for our framework of linear approximation) these func- 
tions will depend on and k'^: — > — A;'^ and (uV) —iuj. 



^We ignore the Weyl anomaly since it is non-linear in the metric perturbations 

*To our understanding, there are no more structures one could possibly add to the expansion (2.23). The only 
tensors with more than four Lorentz indices and linear in h are the ones obtained by applying covariant derivatives to 
the Weyl tensor. 
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The first term generalizes the usual shear viscosity coefficient r/o defined at zero frequency and 
momentum. It also contains the relaxation time term of second order hydrodynamics. The other 
terms (GSFs) are due to metric perturbations, absent in Minkowski space. However, as was pointed 
out in Ref. p9| ], these terms contribute directly to two-point functions of stress tensors, as computed 
from the bulk gravity side. Also the "k" term has been first introduced in Ref. |29| ^. 



From the Minkowski perspective, the physical role of k, p and ^ is not obvious. It is well known 
that the correlators of T^^^ contain not only "thermal" physics but in addition get contaminated by 
the vacuum or zero temperature contributions due to pair production (this is because the underlying 
microscopic theory is a quantum field theory). However, a naive subtraction of T = contributions 
leads to sign alternating results for imaginary parts of the correlators ^] which cannot 

be identified with true thermal spectral functions. This suggests a presence of interference terms 
between "vacuum" and "thermal" amplitudes. 

It is tempting to identify the viscosity term with pure hydrodynamic ("thermal") physics asso- 
ciated with the matter flow, and the GSFs with the non-hydrodynamic or non-matter effects and 
the interference thereof. This conjecture is nicely supported by the ^ term, which at first glance 
does not depend on the fluid's velocity and temperature at all ^. Consequently, when looking at 
the correlators, we will find that in all three channels the contributions due to the ^ term could be 
naturally identified with the vacuum (T = 0) effects. The spectral functions computed from the 
viscosity terms only are positive definite, as they should. 

We would like to comment on the Weyl invariance and non- linear completions. When intro- 
ducing the all order (linearized) hydrodynamics ( 2.23| ) we presented the tensor as transforming 



homogeneously under the Weyl transformation ( 2.16| ). It is obvious, however, that H^'^, as it appears 



in ( 2.23| ), does not have this property. This is because, while the tensor C^°^^^ is Weyl invariant, its 



derivatives are not. Furthermore, higher order derivatives put as arguments of the transport coeffi- 
cient functions also destroy the desired transformation properties. The correct statement is that the 
Weyl invariance is recovered up to non- linear terms, which by themselves are of no interest to us in 
this paper ''. 

It is then a legitimate and interesting question to ask if for any higher order derivative term 
there exists a non-linear completion needed to restore the right transformation property under the 
Weyl transformation. Can it happen that some of the higher order derivatives both in the viscosity 
term and the GSF terms cannot be completed to meet the requirement of the Weyl invariance and 
should be forbidden (similarly to the fate of the bulk viscosity term)? The answer is negative and 
there is no additional selection principle based on the Weyl symmetry. For any number of derivatives 



there exist a non- linear completion with the formal construction given in Ref. [43|. It is based on 
the fact that, instead of the covariant derivative V^, one can introduce an even longer derivative -D^ 
involving the Weyl connection constructed from the field u itself. Any number of these derivatives 
acting on C^J^^ leaves a Weyl-invariant result. This procedure generates non-linear terms, which are 



^In K was introduced as constant. 

®Up to non-linear terms it actually coincides with the stress-energy tensor of the conformal gravity p^ . 
^For the ^ term with constant ^ there exists a well known non-linear completion (see e.g. [^^): under the Weyl 
transformation the tensor Va V/3 C'^°"^^ — 1/2 C^""^^ Ra/i transforms homogeneously. 
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of no interest to us in this paper. For our purposes it is sufficient to know about their existence. 
Note, however, that the procedure of Ref. |43] can be used to reconstruct these non-Hnear terms 

from the higher order hnear terms discussed below. That would certainly provide more insight 
on hydrodynamics at order three and higher. 

If we were not to impose the Weyl invariance, we would introduce another shear viscosity term 
in the expansion ( p. 231) , 

r/2 V" Ua 

This term would normally contribute to the sound channel starting from third order hydrodynamics. 
We would like to argue that this term is in fact forbidden by Weyl invariance. As was explained 
above, in order to comply with Weyl invariance the correct prescription is to use long derivatives -D" 
instead of V". However, the long derivative D° has the property Ua = 0, which eliminates the 
r]2 term. 

The hydro ansatz ( p. 23 ) can be probed by small gravity perturbations. Using linear response the- 
ory we can then compute the retarded correlators in the three symmetry channels (the computation 
is presented in Appendix B). 

• The scalar: 

G^{k,w) = -icori - >,^{w^ + k") - p'-^{w^ - k") + ^^(a;2 - k'f (2.25) 

• The shear: 

G-iK -) = (e + P) ^-^^ - ^---kV2 - pk^k^ - 2.^)/4 + - fe^)/4 ^^^^^^ 

—iuj + r]k'^ 

• The sound: 

G^k, n.) = ie + P) k'-^^n-k--^R.y-2rp.^k^ + -i.^k^ ^^^^^^ 

with 

fj = rj/{e + P); -K = K/{e + P) ; p = p/{e + P) ; C = C/(e + P) ■ (2.28) 

Note that when k = the S0(3) symmetry of the space is restored. Modulo trivial rescaling we do 
indeed observe that the three correlators G"^, G^ , and G^ all coincide: 

2 2 

G^lfc^o = -^G^U^o = -3/4^G^|fc^o - (e + ^ 

At large frequencies w ^ 1, the temperature effects should be negligible and the correlators G 
are expected to coincide with the correlators computed in the vacuum: 

G^{uJ,k)T=o = {e + P) (cj2 - e f ln{k^ - cu^) ; 
G^{uj,k)T=o = -{e + P) k^{uj^ - k^) ln{k^ - lo^) ; 

G^{uj,k)T=o = -{e + P) (4/3) A;2(cj2 _ k'^) in(A;2 - oj^) . (2.29) 
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The asymptotics ( |2.29 ) is indeed observed in the correlators computed from the bulk gravity (see 



previous section). What is interesting to note that the behavior ( 2.29| ) is naturally identified with 
the ^ terms in the correlators, suggesting ^ ~ ln(A;^ — uP') at asymptotically large oj. It is then 
tempting to identify the ^ terms as responsible for the contribution to the correlators of the non- 
hydro pair creation effects, while the k and p terms could be regarded as interference contributions 
between the "vacuum" and "hydro" physics. Within such interpretation it is natural to identify rj 
as purely hydrodynamical effects associated with the matter flow. Thus if one is interested in pure 
thermal/hydrodynamic correlators, one first has to determine r] as functions of momenta and then 
compute the correlators with the GSFs set to zero. 

Despite this nice interpretation of ^ as the pure "vacuum" term, all GSF terms in fact fully 
mix when considered as functions of momenta. If we consider ^ 0, A; ^ oo) asymptotics, all 
correlators tend to behave proportional to fc^ In/c^. From this behavior we can learn about the 
asymptotic behavior of the GSFs themselves 

K^yfcMnA:^, p^VPlnA:^, ^^InA;^. (2.30) 

3. When the bulk meets the boundary: Results 

There should be one to one correspondence between linearized T^^^ and the full set of its correlators. 
Our program is to equate the expressions ( p.25 ,2.26, 2.27| ) for the correlators to the correlators com- 



puted from the bulk gravity. The goal is to invert these equations in order to determine the four 
transport coefficient functions. We have got an apparent problem as we end up having only three 
equations for four unknown functions. This system does not seem to have a unique solution. Despite 
our failure to simultaneously determine all transport coefficient functions, we are able to extract 
them perturbatively in the long-wave limit approximation. 

In the near-longwave limit all of the coefficient functions are expandable in power series ^ 
■q = r/o(l + ir?o,i lo + ??2,o + %,2 + i m,i i^k^ + i %,3 + ^4,0 + ??2,2 k'^ + ??o,4 H ); 

K = Ko{l + iKo,iaJ -I- K2,ofe^ -I- Ko^2W^ + f K2,l W A)^ -|- i Ko,3 "^^ + ■■■)'■, 

p = Pq{1 + ipo^iuj + P2,Qk^ + Pq,2W^ + •••) 

e = Co(i + ^^0,1^ + •••) (3-1) 

Here we explicitly list all terms up to fifth order. The third order coefficients are determined (practi- 
cally all) analytically. The other coefficients are extracted numerically. We achieved a good accuracy 
with the forth order coefficients while the rest have large errors. 

770 = (e + P)/2; r = r/o,i = 2 - ln2; 7?2,o = - 1/2; 

KO = 2r/o; Ko,i = 5/2-2 ln2; po = 4% (3.2) 

The viscosity r/o is of course just ( |1.1D . The coefficient t/q,! is the relaxation time, which within the 



AdS/CFT approach was first addressed in Ref. |44]. It was correctly determined in Ref. Q and 



*We belive this expansion has a finite radius of convergence, The radius of convergence is given by the first singularity, 
which coincides with the first quasinormal mode of the scalar channel. 



- 13 - 



later in [45|. In [^] it was found by looking at the first correction to speed of sound. 770,1 can be 
consistently deduced from any of the three correlators, kq was found also in by matching the /c^ 
term in . Independently and consistently, it can be also found from the shear and sound channels 
(the to k"^ term in the numerator of and the w"^ k"^ term in the numerator of the function G''^). 

The coefficient r/2,0 appears at third order hydro, which was left beyond the scopes of [p9| ]. 
However, this coefficient could be easily read off from the analysis of Ref. [^9|, in particular, from 
the /c^ correction to the diffusive pole in the shear channel. The result is consistent with the k"^ term 
in the numerator of G^ . The coefficient po is deduced from the u! = limit of the function G^ . 
Finally we analytically extracted the coefficient ko,i- This comes from matching the coefficients of 
the u! k"^ in the scalar channel. 

The remaining coefficients were found numerically. Let consider the coefficient ryo,2 as an example 
of our numerical procedure. We were able to get a very accurate fit of the coefficient in front of the 
oj^ term in the expansion of the correlator G^ . This coefficient is then trivially related to 770,2 and 
^0) ^o,ii P01 the latter being all previously determined. The result is 

7?o,2 ^ - 1.379 ± 0.001 - - \ + ^ (3.3) 

where the last expression is our guess for the analytic expression. The error in eq. (^]^), as well as 
other errors quoted below, reflect our confidence in the results provided. 

Despite the fact that we were not able to find a method to extract four unknown coefficient 
functions from three equations, there seems to be a recurrent procedure, which make this task 
possible, at least perturbatively near the long wave limit. The coefficient ^2.0 can be obtained from 
the a; = limit of the sound correlator G^ . Once this one is known, the w = limit of G^ reveals 
the coefficient etc. 

Below we present our numerical results. 
4th order hydro 

r/2,1 = -2.275 ± 0.005; ryo,3 = -0.082 ± 0.003 (3.4) 

5th order hydro 

7/4,0 = 0.565 ± 0.005 ; r?o,4 = 2.9 ± 0.1 ; 7/2,2 = 1-1 ± 0.2 ; (3.5) 

The GSF's coefficients 

K2,o = -1.6 ± 0.05; Ko,2 = 0.04 ± 0.01; 

Ko,3 = -1.95 ± 0.05; ^2,1 = - 1.6 ± 0.2; 

/)o,i = 0.92 ± 0.01; po,2 = -0.68 ± 0.04; p2,o = -0.755 ± 0.005; 

Co = -2.6 ±0.1; Co,i = -1.1 ± 0.2; (3.6) 

To summarize our knowledge of viscosity function r/, we plot it and compare to the IS one (Fig. 
HI). The NS value is, of course, r/ = r/o- For uj, k < 0.4 we can expect up to 15% correction due to 
momenta-dependence of the viscosity function. 
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Figure 3: Viscosity function (divided by rjo): top AdS/CFT; bottom IS. 



On quasinormal modes and analytic structure of the viscosity function 

The quasinormal modes are poles of the retarded correlators. They have been analyzed in all three 



channels in |39]. We would like to argue that entire information about quasinormal modes is coded 
in the viscosity function 77, while the GSF do not have any poles. If this were not true, we would 
observe appearance of identical quasinormal modes in all three channels, which is not the case at 
least for a number of low lying modes. 



*^-)-E^^ (3.7, 



n=0 



We further argue that LOn coincide with the quasinormal modes of the scalar channel (poles of CP"), 
which have been analyzed in the past (see the table below). At /c^ = they can be computed 



quasiclassically for large n (in fact quasiclassics works well down to n = 2) [46| 



u;„(A;2 =0) ~ Wo + n (± 1 - i) (3.8) 

No analytical expression for non-zero k is known. 

The quasinormal modes of the shear and sound channels are obtained from the following disper- 
sion relations. 

- ioj + ri{k^,uj)k^ = 0; - 3uj^ + k^ - Ai'ql.k^ ,uj) uj k^ = (3.9) 
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As well known, these dispersion relations admit hydro dynamic modes as lowest modes in the spec- 
trum. Higher modes will appear as distorted spectrum a;„. Furthermore, the higher the mode the 
less distortion should be present. In other words, the spectra of all three channels will become de- 
generate for high modes. This tendency is clearly observed in the following table copied from Ref. 

m = 1) 



n 


Scalar channel 


Shear channel 


Sound channel 










Re UJn 




1 


±1.954331 


-1.267327 


±1.759116 


-1.291594 


±1.733511 


-1.343008 


2 


±2.880263 


-2.297957 


±2.733081 


-2.330405 


±2.705540 


-2.357062 


3 


±3.836632 


-3.314907 


±3.715933 


-3.345343 


±3.689392 


-3.363863 


4 


±4.807392 


-4.325871 


±4.703643 


-4.353487 


±4.678736 


-4.367981 


5 


±5.786182 


-5.333622 


±5.694472 


-5.358205 


±5.671091 


-5.370784 



Finally we would like to note that from the behavior of the sound dispersion curve [^] one can 
deduce the following asymptotic behavior of the viscosity function 

7?(/c2 ^ ^2 ^ oo) ^ , (3.10) 

2 UJ 

which supports our understanding that it is a falling function at large momenta. 



4. Model for improved causal hydrodynamics 

Though in this paper we do not pursuit any practical applications, we would like to propose an 
improved and causal hydrodynamics for future use by hydro practitioners. 

While we were not able to achieve our prime goal, of deducing the viscosity function in full range 
of frequency and momentum, we were able to get several new coefficients for the small momenta 
expansion. Below we present a resummation scheme similar to IS, which is an ansatz aimed at pro- 
viding a good model for the entire viscosity function. The model is constructed with the requirement 
of causality built in. 

Causality implies that the imaginary part of the poles is always negative and the function vanishes 
at infinite frequencies. This is equivalent to the validity of the dispersion relation: 

J 111 I uj' — OJ 

In addition, in order to relate the viscosity function to the thermal spectral functions, we require 
that both real and imaginary parts of it remain positive for all values of momentum and frequency. 

Similarly to the IS model, we take a Pade-like resummation ansatz which reproduces all low 
momentum coefficients in the expansion. 

3 ^ 

ryweP =^oE ,^ + 5.fc2_,^ (4.2) 

i=l 
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This ansatz has three pure imaginary poles and it reproduces exactly eight first coefficients in the 
expansion (|3.1|). 



di = 0.736, ai = 0.72731, 61 = 0.3263 ^2 = 2.1, 02 = 0.10618, 62 = 0.3042, 
^3 = -2.1016, 03 = 0.10620, 63 = 0.3038. 
The resummed viscosity function is plotted in Fig. ^. This model could be further improved by 




Im [77] 




Figure 4: Viscosity function (divided by r]a): the model 



accounting for the asymptotic behavior ( 3.1C| ) as well as for information about quasinormal modes of 
the scalar channel. The second and third poles practically cancel each other. Despite the fact that it 
does not accurately reproduce the expansion, it turns out to be a very good approximation to retain 
only one pole, similarly to IS but with three-momentum dependence. 

Vo 



Vmodel 2 



(4.3) 



Within about 10% accuracy (and in some regions with much better one) the second model is equiva- 
lent to the first one. Since the entire effect of momenta-dependence is not expected to be very large, 
the second model should be more than sufficient for any phenomenological applications. We note 
that the group velocity for the sound mode computed within this model is always smaller than one, 
confirming causality of the model. 

The viscosity function can be Fourier transformed into the memory function 

-^^t + ikx ^^j^2^^-^ (4.4) 



D{x,t) = J diodrke 
which leads to the following expression for the dissipation tensor H: 

= - 2 dt' J d^x' D{x -x',t- t') u^ix', t') 



(4.5) 



Performing the Fourier transform explicitly we obtain 
We remind the reader that r/2,0 is negative. 



1 Vo 
2 V2 Vo,i 



-Vo,i 



3/2 



(4.6) 
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5. Summary and Discussion 



In this paper we initiated a study of all order velocity gradient expansion of linearized relativistic 
hydrodynamics near equilibrium. The research was carried out within the = 4 SYM theory at large 
Nc- More specifically, we parameterized the energy- momentum tensor of the theory in terms of four 
momenta-dependent functions. These functions generalize the notion of the usual constant transport 
coefficients, such as viscosity, into momenta dependent ones. We then attempted to determine all 
four functions based on the information on retarded correlators of two stress tensors. The latter 
were computed via the AdS/CFT prescription for computing retarded correlators from bulk gravity 
waves. 

Out of four transport coefficient functions, which we introduced in ( |2.23| ), r/ appears as a coef- 
ficient of the operator constructed from velocity gradients and is a generalization of shear viscosity. 
The remaining three coefficients (GSFs) arise as coefficients of 4d metric perturbations and appear 
in front of three operators involving the curvature Weyl tensor. 

In this paper we extended the previous knowledge of hydrodynamic transport coefficients at first 
and second order to some higher order coefficients. We were able to find only those, which contribute 
to linearized hydrodynamics. We gave analytic values for two coefficients of the third order hydro. 
We provided very accurate numerical estimates for two coefficients of the forth order and one of the 
fifth. In addition, we introduced and determined several new coefficients associated with the GSFs. 

To illustrate the effect of the higher order terms in the viscosity function, we compute the sound 
dispersion curve by solving (perturbatively) eq. (|3.9| ) to the order k^: 



(l + Q - - 0.088/^4^ -i^ _ ^(4 - 81n2 + In^ 2) - Q.lbk"^ 

(5.1) 



As we have already emphasized in Ref. |31], the sound width gets negative corrections from higher 
order terms. This is in sharp contrast to the IS model, which leads to a qualitatively opposite effect, 
with the correction being positive 

Based on the new information about higher order terms in the expansion of the viscosity function, 
we have proposed an improved causal IS-like (single pole) hydrodynamics, which we hope can be 
used by hydro practitioners. Compared to IS, this model emphasizes the importance of the space 
momentum dependence of the viscosity function. It leads to qualitatively different predictions, as 
seen from the width of the sound pole. On the basis of this example, we cautiously suggest that the 
results based on the IS theory might be in fact less reliable than it was previously thought. We also 
propose to exploit our improved model for non-linear phenomena (such as Bjorken expansion and 
elliptic flow) even though such applications have no theoretical justification. 

Admittedly, the problem we had set up was not yet fully solved by the present paper. Using a 
perturbative procedure, we found several new higher order (constant) transport coefficients, either 
analytically or numerically. Nevertheless, a generic problem remains to be solved: the correlators we 
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used as our input seemed not to be sufficient to determine the four transport coefficient functions 
parameterizing all relevant kinematic structures. While it is possible to follow the iterative approach 
used by us to determine more coefficients in the expansion of the transport functions, we have no 
proof that this procedure will actually work at all higher orders with unique results. It is not excluded 
that some additional inputs (apart of the correlators) are required in order to solve the problem in 
full. 

As an alternative approach to the problem, one may switch from solving the bulk equations for 
gravity waves - the basis for computing the correlators - to a membrane paradigm-type approach 
based on vibrations and translations of the horizon, as it is done in Ref. and the follow on papers 
1 47 , p4| . This approach provides a quite general procedure to derive the next order derivative terms for 



boundary hydrodynamics. In this approach, the boundary energy-momentum tensor with appropriate 
gradient corrections is obtained through the usual holographic renormalization procedure, with the 
bulk solution reflecting perturbation of the near-horizon "membrane" . If the boundary metric is not 
taken as flat Minkowski (as it was done in |^0|), but rather as a slightly perturbed one, the method 
of [^] would reveal the GSFs alongside the viscosity function. Furthermore, the approach of [^] 
has a potential to determine not only linear but also non-linear terms, the latter being beyond the 
scope of our present paper. Here we obviously mean third and higher order hydrodynamics. We have 
not pursued this direction, but believe it is worth studying it as it is important to learn about the 
gradient structure as a way to understand the non-equilibrium effects in plasma. 

An important general problem is a separation between the hydrodynamic (thermal) physics 
associated with the matter flow and the vacuum (zero temperature) effects associated with the 
pair production, as both contribute to the retarded correlators. We hope that we proposed the right 
approach to it, by identifying the different roles played by the viscosity function and the GSFs. While 
the former is purely "hydrodynamical" , the latter includes non-thermal physics and interference. This 
separation of roles is very plausible, supported by the results at hand, but it was not proven by us 
in general. We have argued in the text that the pole structure of the correlators is entirely included 
in the viscosity function, while the GSFs have no poles. The overall role of the GSFs is somewhat 
unclear. On the one hand, they are formally introduced in ( 2.23| ) as a response of the fluid to 



external gravitational shakings. On the other hand, from the analysis of the correlators we identify 
the GSFs as being responsible for flat space non-hydrodynamic (non-thermal) effects associated with 
pair production of the underlying microscopic field theory. The metric perturbations in effect mimic 
non-hydro physics. 

Are the GSFs relevant for RHIC experiment? We believe the answer is "no" . We had to deal with 
them only because the correlators used for our analysis contain both thermal physics and vacuum 
effects (such as pair creation). One, of course, could propose another type of experiment, in which 
plasma would be exposed to a real gravitational wave. In this type of experiment, the GSFs would 
determine the physical response of the fluid. 

The so-called contact (or Schwinger) terms are QFT phenomena originating in UV. One could 
suspect that the GSFs originate from those. However, if this were the case, the only effect they 
would produce is to shift the correlators by finite order polynomials in momenta. This is not the 
case, however. From the explicit expressions for the correlators ( |2.26| , |2.27| ) one can see that the 
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numerator terms involving the GSFs cannot cancel the corresponding poles. Thus they include more 
physics than just the contact terms part of which is not coming from the UV. 

Last but not least, it remains to be seen how relevant the effects of momenta-dependent viscosity 
are for realistic applications of relativistic hydrodynamics to heavy ion collisions. Our previous paper 
1 31] argues that they might be quite substantial at the early times of the collision. As we explained in 
the Introduction, recently the phenomenology shifted to the "fate of small initial state fluctuations" , 
related with the conical structure and "ridges" . Although we have not applied our results, it is clear 
that such flows would be sensitive to higher gradients, as the size of those fluctuations is an order of 
magnitude smaller than the nuclear size associated with radial and elliptic flows studied before. In 
agreement with our proposal |31|, the viscosity function (its real part) is a decreasing function both 
of frequency and momenta. This behavior might be the reason behind the low viscosity observed at 
RHIC. It may also explain the exceptionally good survival of various hydrodynamic flows, particularly 
the sound waves. 
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A. Appendix: Analytic expansions of correlators 

In this Appendix we present analytic expansions of retarded correlators at small frequency and 
momentum, as computed from the bulk gravity. The expansions are obtained following Appendix of 



Ref. ]|29[] where a perturbative approach to solving eq. ( p.2| ) is set. We reproduce and extend their 



results to include some of higher order terms. 
Scalar channel 

ln2 , /3 1n2 \ ^ , In^ 2 , . In^ 2 , , /5 ln^2\ . 

A = 1 + i a; + ln2 l] uj^ + ln2k^ - i ujk^ + k^ 

2 V8/ 2 2 14 6/ 
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1,9 1 /3 ln2\ ,4 /I ln2\ , ln2 ln2 , , . , « , 

S = -A;2 + i-u^ + _ (^^ " — j ^ " ^ — ^^ + — (3-ln2) /c^ ... (A.l) 



For the retarded correlator we obtain 



— ^— = = --A:2-iiw--(ln2-l)cj2-i (3 - 41n2) A:^ + Hn2LjA;2-ln22A;^ 

(e + P) A 2 2 2 ^ ' 4 ^ ^ 



(A.2) 



Shear channel 



1 



1,4 ln2 



In 2 



A = uj + i — k + i — k + w k'^ + i lj" 



B 



2 4 

.2\ 



« 2 . 2 — In 2 



■ UJ 



.k\..) 



(A.3) 



The correlator reads 



1 



G 



D 



B i k'^/2 [1 + i (2 - In 2) w - k'^/2 . . .] + lo k'^/2 + 



k^ A 



CO + iA;2/2[l + i(2 - ln2)6j - ^2/2 + 



(A.4) 



Sound channel 



^ = 8 



In 2 In 2 ' 

9uj^ - 3k^ + iGujk^ + (ln2 - 4) A:^ + 9 ln2(^ 1)uj^ + 3 ln2(2 ^)^^'^^ 



5 = 3 [I8u;2 - SOk"^ + a2Lj/c2 + 2(5 ln2 - 12)k'^ + 3 ln2(12 - 5 ln2)u;^A;2 + 

+ 9 In2(ln2 - 2)w^ ...] (A.5) 



The analytically controlled part of the sound correlator 
1 „c /4 k^ 



G- 



B 1 29A:^ - SOk'^LV^ + 90;^ 3 
3(^2 - A;2)2 1 + 12 (/fc2 _ a;2)2 4 

/c^ + i2[l - iu;(ln2 - 2) + •••]a;A;2 _^ 2a;2A:2 ... 
3lj2 - /c2 _^ i2a;A;2 [1 - iu;(ln2 - 2) + • • •] 



+ 1 



(A.6) 



Here we used eq. ( p. 10 ) 



B. Appendix: Correlators from generalized hydrodynamics 



In this Appendix we compute the retarded correlators from the hydrodynamic ansatz (|2.23| ) using 4d 
metric perturbations. The non-perturbed space has the Minkowski metric g^'^ = diag{—l, 1, 1, 1}. 
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Scalar channel 

The perturbation is h = hxy{z, t). The fluid remains at rest. We first compute Christoffels coefficients 
and Riemann tensor 

r* = = = i /) ■ = vy = — = ^h' 

xy ty tx 3 y 2 

Tjt _ TJX — jyy — }_ i . — jyy = _ = — h" ■ 

^xty -'^tty ^ttx 2 ' -^zzy ^zzx ^xzy 2 ' 

R* _ nx = = — Ti^ = T}^ = = }l h' CP, ^\ 

^xzy -'^tzy ^tzx ^xzy ^zty ^ztx K'^-'-J 

The only non-zero component of the Ricci tensor is Rxy while the scalar curvature is zero 

Rxy = \ih- h") R = (B.2) 

The relevant non-zero components of the Weyl tensor are 

Ctxty — Cfxyt — Cxtyt — Cytxt — Cxzyz — Cyzxz — ^ ~l~ ^ ) 

Cxzyt — Cyzxt — Cxtyz — Cytxz — 2 ('^'^) 

The xy component of the stress tensor reads 

(T^y) = -Ph -r/h- ^K[h + h"] + p^ih -h"] - C^ih -2h" + h""] (B.4) 



In momentum space this becomes 



(rji<xy\^ 



h{k,w) (B.5) 



From Eq. (|B.5| ) one can read off the correlator = G^y^y 



G''y'-'y{k,w) = P - iuT] - K^iw'^ + k^) - p'-^{w^ - k^) + (^{uj^ - k^'f (B.6) 

The retarded correlator G^y^y = G^y^y — P. The transverse static susceptibility is momenta 
dependent and is given by the functions k and ^: 

X^{k) = - Kik, 0) k^/2 + 0) k^i (B.7) 

Shear channel 

The perturbation /i = /ij^(2;, t). The fluid's four velocities is = (1, u, 0, 0) and = {—l,v+h,0,0) 
The Christoffels coefficients and the Riemann tensor are 



■pa; 

'-tz — ~ 


r* = - 

^ xz 


^ xt - 


2 ' 


J- tt 


^txtz — 


nt _ 
^xzt 


ryx 

^ttz 


— L)^ — 

- - ^tzt - 


~ ^txt 


TJZ _ 

^xzt 


~ ^xtz — 


nx 
^ztz 


= ~ ^txz = 


-2^ 



Ik' 

2 



(B.8) 
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The non-zero Ricci components and curvature are 



Rxz = ~2^'' ~ ~2^"' R = (B.9) 



The relevant non-zero Weyl components read 



Cxtzt — Cztxt ~ ^ ^' ' Cxzzt — Cztxz — -^h" (B.IO) 

The components of the stress tensor 

(r«) = e; T*^ = (e + P)v + Ph 

(T^-) = -r^v' + K^h' + (h'" - 2h') + C^ih - h'") (B.ll) 
Equations of motion relate the metric perturbation h to the induced three- velocity v: 
dt{T'-) = {e + P){v + h); 

dz (T-) = -rjv" + K^h" + ih"" -2h" + (h - h"") (B.12) 

which leads to the relation (in momentum space) 

, iu - iRuk'^/2 - pe{e - 2a;2)/4 + i^uk'^iuj'^ - k'^)/4 ^ , 

V = h : _ 2 (B.13j 

Substituting this relation back into the expression for (T*^) we can read off the correlator = G*^*^ 
Gt.t. = u + P) ^fe^ - ^^^fcV2 - -pk\k^ - 2u;^)/4 + i^^jk^ju^ - k^)/4 _ ^ 

Note the appearance of the extra terms proportional to the GSFs in the numerator, whereas in the 
normal diffusion scenario the residue is usually given by the viscosity only. 

The retarded correlator G*^*^' = (5*^*^' + e. The shear static susceptibility x^'- 



X^ik) = (e + P) 



1 _ -P{k,0) ^2 
4f?(/c,0) 



(B.15) 



Sound channel 



The perturbation which generates sound is h = htz{z, t). The fluid's four velocities is u^^ = (1, 0, 0, v) 
and -u^ = (—1,0,0,^-1-/1). 

Christoffels and Riemann are 

r'L = ~ h' ; = h; P-tztz — ~ P-ttz = P-zzt — ~ P^ztz (B.16) 

Contrary to the cases of tensor and shear perturbations, the sound perturbation has a nonvanishing 
scalar curvature. 

Rzz = -Rtt = -h'; R = - 2h'; C,tzt = \ h' (B.17) 
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The relevant components of the stress tensor 



(r«) = e; {T'') = {e + P)v + Ph; 

(T-) =P- n^v' + - p^h' + h (B.18) 



Equations of motion can be solved for v relating it to the perturbation h 



k"^ — 3u;2 — Aifjuk"^ 



Substituting v back into the expression for T*^ we can read off the the correlator = G*^*^ 

^ — 2Ku^k^ — 2ipu' 
k"^ — 3a;^ — Aifjujk^ 



G = {e + P) -2 —2 jT^-To ^ l^-^") 



The retarded correlator G*^*^ = Qtztz _|_ ^ 
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